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We calculate frequencies and damping rates of the lowest collective modes of a dilute Bose 
gas confined in an anisotropic trapping potential above the Bose-Einstein transition temperature. 
From the Boltzmann equation with a simplified collision integral we derive a general dispersion 
relation that interpolates between the coUisionless and hydrodynamic regimes. In the case of axially 
symmetric traps we obtain explicit expressions for the frequencies and damping rates of the lowest 
modes in terms of a phenomenological collision time. Our results are compared with microscopic 
calculations and experiments. 

PACS numbers: 03.75.Fi,03.65.Db,05.30.Jp,32.80.Pj 



I. INTRODUCTION 

Above the Bose-Einstein transition temperature both coUisionless and hydrodynamic modes can exist in a trapped 
gas. In the coUisionless regime the frequency of the mode is large compared with the interatomic collision frequency, 
and the mean free path is larger than the wavelength of the mode and the dimension of the cloud. In the hydrodynamic 
regime the opposite is true. Previous theoretical studies have focused on the coUisionless and hydrodynamic 
regimes Analysis of recent measurements of frequencies and damping rates of collective modes in trapped 

Bose gases above the Bose-Einstein transition temperature P-pd| indicates that the experiments were performed 
under conditions intermediate between the coUisionless and hydrodynamic regimes and in Ref. [Q a simple 

interpolation formula for the intermediate regime was proposed. The aim of this paper is to calculate frequencies and 
damping rates of modes in the intermediate regime, and to examine the validity of the interpolation formula. 

The difficulty in this problem arises from the inhomogeneity of the system. Unlike the case of homogeneous systems 
conditions may be hydrodynamic in the center of the cloud and coUisionless near its surface. In a trapped gas the 
coUisionless limit is always well-defined in the sense that the frequency may exceed the collision rate everywhere, 
but hydrodynamic conditions can only be achieved in a limited region of space, not in the outer parts of the cloud. 
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Therefore one would expect a single relaxation-time approximation to work less well than for a homogeneous system. 

At temperatures above the transition temperature the gas is dilute in the sense that the mean interatomic spacing 
is much larger than the interatomic scattering length. Mean-field effects can be neglected since the potential energy 
due to interactions, nUo, is much less than the thermal energy ksT. Here n is the number density, Uq is the effective 
two-body interaction (which is given in terms of the scattering length a as Uq — 47raft^/m, where m is the atomic 
mass), kg is the Boltzmann constant, and T is the temperature. Then the atoms can be described by a distribution 
function that obeys the Boltzmann equation [Q . We solve the Boltzmann equation with a trial solution appropriate 
to the collective modes being studied, using a relaxation time approximation to the collision integral. 

In the next section we derive a general relation for mode frequencies that interpolates between the collisionless and 



hydrodynamic regimes in an anisotropic trap. In section III we consider the experimentally relevant case of axially 



symmetric traps and calculate frequencies and damping rates of the lowest modes. In section IV we discuss the validity 
of our approach by comparing our results with microscopic calculations and with experiment. Section ^ contains a 
summary of our conclusions. 

II. THE BOLTZMANN EQUATION AND COLLECTIVE MODES 

In this section we derive the dispersion relation of the lowest modes of a Bose gas confined by an anisotropic 
harmonic potential given by 

V{r) = ^m{u;lx' + cj^y' + u^jz^), (1) 

where t^i , , and W3 are the characteristic frequencies of the trap in the three directions. 

The dynamics of the dilute Bose gas above the transition temperature can be described by a semi-classical dis- 
tribution function /(p, r, i), where p is the particle momentum. The distribution function satisfies the Boltzmann 
equation 

9/1 fdf 



^ + -p.V./-VrV^.Vp/=^) . (2) 



coll. 



Here {df/dt)^^/^^ is the contribution of collisions to the rate of change of /. For two-body scattering it is given by 

(^^^%^) =-(^ / ^Pi / d^|v-vi|[/(p,r)/(pi,r)(l + /(p;r))(l + /(pi,r)) 

-/(P;r) /(pi,r)(l-f/(p,r))(l + /(pi,r))]. (3) 
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The collision integral involves scattering processes with two incoming particles of momenta p and pi (and velocities 
V and vi ) and two outgoing particles of momenta p' and pi , respectively. The differential cross section is denoted by 
da. 

We shall linearize the Boltzmann equation in small deviations of the distribution function from the equilibrium one. 
This can be achieved simply by inserting the following expression 

/(p, r) - (p, r) + (p, r)(l + /(") (p, r))^(p, r) (4) 

in Eqs. ^ and (||), where V'(Pjr) describes the deviation from equilibrium, and /(°) is the equilibrium distribution 
function 

/(")(p,r) = {exp[(pV2m + V^(r)-M)/fcsT] (5) 
To first order in -0, the Boltzmann equation assumes the following form 

^ + ^p • Vr - VrF • Vp^ V(P, r) = -imp, r)], (6) 
where / is a linear operator given by 



(2^ / / ^'^l^-^il /^°^(Pi,r)(l + /("Hp;r))(l + /(")(p'i,r))^^^^i^ 



r 



(7) 



[ij{p, r) + V^Pi, r) - ?A(p; r) - ?A(pi, r)] . (8) 

Since the scattering processes considered here conserve the number of particles, momentum, and energy p^ , the 
zeroth, first, and second moments of the collision integral, calculated by multiplying it by 1, p, and p^, respectively, 
and then integrating over p, must vanish. The quantities 1, p, and are called the collision invariants. 



A. Dispersion Relation for Low-Lying Modes 

The lowest collective mode corresponds to the center of mass motion (the dipole mode). The frequency of this 
mode equals cui {i — 1, 2, 3) and does not depend on interatomic collisions. Therefore this mode will not be considered 
further in this paper. We shall instead focus on the higher modes which have a frequency 2uji in the coUisionless limit. 
In the hydrodynamic limit one can show that the fluid velocity u associated with these modes is given generally by 
u = {ax, by, cz), where a, b, and c are constants 1^-0]. This expression is general in the sense that it does not depend 
on temperature or the anisotropy of the trap. In spherical traps these modes correspond to the monopole mode (or 
breathing mode) characterized by a = 6 = c, and the quadrupolc modes for which a + 5 + c = 0. 
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Let us attempt to find a solution ip that corresponds to the lowest modes of the system described above. The 
deviation function ■0 corresponding to a flow with drift velocity u is proportional to u • p. Acting on u • p with the 
left hand side of Eq. (^) will thus introduce terms like x^, pi etc., and therefore we adopt as our trial function the 
general form 

-0 = e"'"*[aia;^ + bixp^, + cipl + 02?/^ + b2ypy + C2pl + 032;^ + 632^^ + cs.pl]. (9) 

Here the terms x^, y^, and correspond to changes in the number of particles, and xp^, ypy, and zpz correspond 
to changes in the local momentum density, and hence all these terms are collision invariants. On the other hand the 
terms p^., Py, and p^. are not collision invariants since only their sum p^ = p^. + Py + pI, which is proportional to the 
kinetic energy, is. Mode damping therefore arises from only these latter terms. 

In the relaxation time approximation the collision operator is associated with a mean relaxation time T(r), which 
depends on the position r. One is restricted, however, in that for any approximate expression for the collision 
operator /, the above-mentioned conservation laws should be satisfied. This can be ensured by requiring the collision 
invariants to be eigenfunctions of the collision operator with eigenvalue zero. However, for terms such as p^, which are 
not collision invariants and in momentum space have the symmetry of 5^2™! where Y"^ are the spherical harmonics, 
we have 

i\pl]^^^^ipl-pVs), (10) 

and similarly for Py and p^.. The term p^/3 in Eq. (|l^) ensures that I[pl + Py + pI] = 0. 

Since the scattering rate depends on the local number density p^ , the mode damping involves a spatial average 
of the relaxation time. (The simplest functional dependence one can assume for the collision rate is that of the 
equilibrium distribution fmiction.) To obtain collective modes, we insert the trial function (|^) in the Boltzmann 
equation and then take moments of this equation by multiplying it by /'■''-'(I + f^^^) and by each of the terms in 
-0 and then integrating over p and r taking into account Eq. ([l^). For each moment we get a separate equation and 
thus obtain nine coupled equations. The condition for the existence of nontrivial solutions results in the following 
dispersion relation for the frequency of the collective modes: 

+ li (c^6 _ 10cj2^4 ^ 32^4^2 _ 9g^6^ (11) 

where 
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2 2,2,2 



4 2 2 2 2 2 2 

= ^1^2 + ^l^Z + ^2^3: 

and 
with 

l_ /drr-i(r)/dp/W(l + /W)£^ 

r /dr/dp/(")(l + /(o))e2 ' ^ > 

where e — /2m is the kinetic energy. Equation ( ITTI ) is our main resuh, which describes the frequency and attenuation 
of collective modes in the coUisionless and hydrodynamic regimes, as well as in the intermediate regime. 



III. MODE FREQUENCIES AND DAMPING RATES FOR A TRAP WITH AN AXIS OF SYMMETRY 

In traps with an axis of symmetry the oscillations in the axial and radial directions are decoupled and therefore the 
projection of angular momentum m about the symmetry axis of the trap is a good quantum number and can be used 



to label modes. In this case the modes described for a general trap in subsection II A have the following three forms: 
i) The oscillations in the radial and axial directions are in phase (In the spherically symmetric case this corresponds 
to the breathing mode, a = b = c). ii) The oscillations in the two radial directions are out-of-phase with those in 
the axial direction { a = b = — c/2 for the case of spherical symmetry), iii) The oscillations in the two perpendicular 
radial directions are out-of-phase with respect to each other and there is no oscillation in the axial direction (a = — &, 
c = 0). Both of the first and second modes have ni = 0, while for the last one m — 2. For the rest of this paper we 
denote the first mode as the 0+-mode, the second one as the 0^-mode, and the third one as the 2-mode. 



A. Frequencies and Damping Rates 

The roots of Eq. ( p^ ) can be obtained analytically as functions of r and the characteristic frequencies uJa, Wfj, and 
ujc, but they have a complicated functional dependence on the parameters. However, for experimental traps with 
axial symmetry characterized by uji = uj2 = and uj^ — Xloq, where A is the anisotropy ratio, the roots simplify 
considerably. In this case the dispersion relation ( pi] ) takes the form 

LU [{uj^ - 2ujI) - iujT - 4cj^)] [{lu^ - r_wg) (w^ _ r+wg) - iivT {uj^ - 4w^) (w^ - 4A^cj^)] = 0, (13) 
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where 

r±= ^ + 4A^ ± ^ V25 - 32A2 + 16A4. (14) 

Note that the dispersion relation in this case factorizes into two terms (apart from the lj factor), one of which is 
independent of A. This factorization can be understood by finding the eigenmodes that correspond to the roots of 
(pi]). As we shall discuss in more detail below, the A- independent term corresponds to the 2-mode, in which there 
is no oscillation in the z-direction. The other two modes, corresponding to the A-dependent factor, are the 0"*"- and 
0"-modes. There are three modes of frequency w = in the coUisionless limit which correspond to thermal expansion. 
The first corresponds to constant temperature change in all directions. This mode is associated with the over-all 
factor u) in Eq. ([l3|). The second mode has m — and frequency uj — —ijlT in the coUisionless limit (u^t^T S> 1), and 
Lo — —ijT in the hydrodynamic limit (wqt <C 1). The third mode has m — 2 and the same values of frequency in the 
coUisionless and hydrodynamic limits as for the previous mode. 

We start the analysis of (^3|) by looking at mode frequencies in the coUisionless and hydrodynamic limits. In the 
coUisionless limit, wr ^ 1, we get: 

a; = tJc = 2tJo, 2ijJo, 2A(jJo, (15) 

while in the hydrodynamic limit, lot <C 1: 

UJ = UJH^ \/2tJo, ^/rI^^n, \AT^o, (16) 

where wc denotes the frequency in the coUisionless limit and loh is the corresponding frequency in the hydrodynamic 
limit. Thus we obtain the simple picture that the three coUisionless modes of frequencies 2a;o, 2ajo and 2Awo correspond 
to the three hydrodynamic modes of frequencies v^'^Oi ^/t^^o and ^Jr^uo^^ respectively. (This is true only for A > 1. 
When A < 1 the mode of frequency 2Awo corresponds to the mode of frequency y^rTcjo instead of that of frequency 
y/r^ujQ.). We note here that our result ( p^ ) agrees with previous works based on either the conservation laws or 
taking moments of the kinetic equation (H . 

To obtain frequencies of modes and their attenuation in the intermediate regime, we insert oj = ujr — iuji in Eq. 
(|l3|), and thus obtain two coupled equations for the real and imaginary parts. In Figs. 1 and 2 we plot the frequency 
and the damping rate of the 0^-mode for values of A that satisfy the experimentally relevant condition A <C 1. 

It is instructive, also, to investigate the functional dependence of the damping rate on A and t in the coUisionless 
and hydrodynamic limits. This behavior can be deduced using the property that in both of these limits the imaginary 



part uji is small, and thus one can expand Eq. ( [T^ ) in uji. Let us first consider the dispersion relation corresponding to 
the first factor in this equation. In the coUisionless limit, lot ^ 1, we can expand this factor to first order in 1/{ujt), 
and write the solution as 

w = 2^0 - ir^, (17) 
where is the small damping rate associated with the 2-mode. To first order in 1 /r we obtain 

n = lr-\ (18) 



Similarly, for the other two modes associated with the A-dependent factor of (^^ we obtain 

1 
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and 



= T^r-\ (19) 



n-=lr-'- (20) 



In the hydrodynamic limit, ujt ^ 1, we expand in lut and then insert lu = V^ujq — to get 



together with 



r'^ = ^or, (21) 



, (r+-4V)(r_-4) 

^ 2(r_-r+) (2^) 



and 

, _ (r_-4A2)(r_-4) 
^0" - 2(r+-r_) 

We note that the damping rates in the coUisionless limit are proportional to 1/t and independent of A. 

Finally, we end this section by discussing some interesting features of the hydrodynamic damping rates Fp+ and Fg_ 
given by Eqs. ( p2|) and (|2^). This can be done by plotting these rates as a function of A as shown in Fig. 3. From this 
figure, we first note that F^^ vanishes for A = 1. This is in agreement with the well established fact that in an isotropic 
harmonic potential the breathing mode does not relax, as was first shown by Boltzmann |jl^. We also note that Fq. 
has the maximum value I.ISwqT for A ~ 1.20. In the limit A — > (cigar-shaped cloud; a quasi one-dimensional 
system) the damping rate vanishes, while in the limit A — > oo (disk-shaped cloud; a quasi two-dimensional system) 
the damping rate approaches a constant value equal to {3/A)luqT. 
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IV. DISCUSSION 



The damping rates found above were all expressed in terms of the phenomenological time r, Eq. ([T^). We shall 
now identify the damping rates found above with the precise expressions obtained by a microscopic calculation. If we 
assume that the collision rate is proportional to the local density and take the classical limit, ^ 1, the spatial 
average entering Eq. ( [T^ ) may readily be performed, resulting in r = 2V2t(0), where t(0) is the relaxation time 
at the center of the trap. The damping rate, ro+, for the O^-mode, which is given by Eq. (|l9|), therefore becomes 
l/(24\/2r(0)). In Ref. § this damping rate was calculated from a variational solution to the Boltzmann equation in 
the collisionless limit, with the result 

where the quantity in the parenthesis results from including degeneracy effects to lowest non-vanishing order. Here 
TvariO) is the variational viscous relaxation time given by t~J^{0) — {Ay/2/5)n{0)av where n(0) is the central 
density, a = Sna^ is the total cross section, v = y^SksT/ {nm) is the mean thermal velocity, and Tbec is the Bose- 
Einstcin transition temperature. The use of an improved trial function changes the number 1/(24V2t(0)) in Eq. ( [2^ ) 
by only ~ 1% A comparison between the damping rate Eq. ( p^ ) and Eq. (^9|) thus reveals that 

r(0) = (^1 + ^C(3) (^)') W(0). (25) 

The same relationship, Eq. (|2^), applies to the other two modes as well. For the two experiments of Refs. [||J|] our 
phenomenological time, t(0), in the center of the cloud is therefore given by t(0) ~ 0.97Ti,ar(0) for the experiment 
performed at T ~ 2Tbec, and t(0) ~ 0.82Ti,ar(0) at T ~ Tbec- 

In the hydrodynamic limit, an accurate microscopic calculation of the collision rate is difficult to obtain due to 
the failure of the hydrodynamic conditions near the surface of the cloud. This is because the mean free path is too 
large for hydrodynamics to be valid. To solve this problem a cut-off procedure was introduced in Ref. [Q, taking the 
hydrodynamic description to be valid to a distance from the center of the cloud such that an atom incident from 
outside the cloud has a probability less than 1/e of not colliding with another atom. This cut-off procedure led 
to a logarithmic dependence of the damping on the quantity a — an{Q)^JkBT / (2mAwQ). The leading term of this 
logarithmic dependence gives a damping rate proportional to [Ina]^/^ [Q. The trial function (||) is not capable of 
describing accurately the hydrodynamic limit, since the very concept of hydrodynamic flow loses its meaning in the 
surface region. Expressed in other terms, the deviation function has a very different form in the outermost part of 
the cloud than it docs in the center, and this effect is not captured by our simple trial function. 
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The microscopic calculations performed in Refs. accounted for the damping rates only in the coUisionless and 
hydrodynamic limits. To compare with experiments, which are performed under intermediate conditions, a simple 
phenomenological interpolation formula was suggested Q. In the following we discuss the validity of this formula. In 
particular we show that the interpolation formula is identical to our result for the 2-mode. Our results for the 0+- and 
0~-niodes cannot in general be expressed in this simple form, but show a close resemblance, as we shall demonstrate 
below. The interpolation formula suggested in Ref. [Q is 

^2^^2^i^i^^. (26) 

1 — iU!T 

By construction, the limits lut 3> 1 and lut <^ 1 yield the appropriate frequencies in the coUisionless and hydrodynamic 
limits, respectively. The 2-mode is A-independent, thus one can easily put the corresponding dispersion relation 



(resulting from equating the cubic factor in uj in Eq. (13) to zero) in the same form as (p6[). For the -mode we 



consider the two limiting cases A <C 1 and A 3> 1, and the corresponding dispersion relation can be written as 

1 — 2Cj(6r/5) 

- + = 0, A » 1). (28) 

1 — za;(4r/3) 

The formula corresponding to the O^-mode in these limits takes the same form as (^), but differs slightly from it, 
when A is comparable to 1. To illustrate how well the simple interpolation formula, Eq. (^6|), agrees with our results 
obtained from the solution of Eq. ( p^ ) we compare in Fig. 4 our calculated damping rate of the 0^-mode, for the 
case A = VS, with the interpolation formula. 

From experiments one obtains independent measurements of the mode frequency uj,, and its damping rate Wj. Since 
our calculated values of these quantities are functions of the collision rate, we can eliminate r and obtain a relation 
between the real and imaginary parts. This allows us to calculate the damping rate of a mode without knowing about 
the characteristics of the collision processes leading to damping. Substituting uj = Ur — iuji in Eq. (^) and eliminating 
T from the two resulting equations leads to 



4 



2 



(29) 



It is interesting to note that this expression is also what one obtains from the simple interpolation formula (p9). In 



fact the expression Eq. (E9h is valid for any number replacing the ratio 6/5 in Eq. (27). In Fig. 4 we use Eq. (|29|) to 



plot LOi as a function of Ur- Two values of A are used, namely 0.076 and 0.074, corresponding to the two experiments 
of Refs. It can be clearly seen in this figure that the conditions of these experiments are intermediate between 

the collisionless and hydrodynamic limits. The deviation of Eq. ( |29| ) from the exact solution obtained from Eq. ( p^ ) 
is shown in Fig. 5 where we plot uji as a function of ujr for A = \/8- 

The present method of obtaining the frequency and attenuation of collective modes can be extended to higher 
modes of frequency luji, where I is a positive integer, by employing trial functions that differ from the one we have 
used in Eq. (|9[), for the 2tJi-modes. 



V. CONCLUSION 



Starting from the semiclassical Boltzmann equation and using the relaxation time approximation we have derived 
expressions for the frequencies and damping rates of the low-lying oscillatory modes of a dilute Bose gas above 
the Bose-Einstein transition temperature. These expressions give frequencies and damping rates as functions of a 
phenomenological collision time and hence describe modes in the regime intermediate between the collisionless and 
hydrodynamic regimes. Our results justify the use of a simple interpolation formula for connecting the hydrodynamic 
with the collisionless regime. The approach used in this work may be readily generalized to higher-frequency modes. 
After the completion of the present work we received a preprint by Guery-Odelin et al. who consider the O"''- and 
0^-modes in the classical limit and obtain conclusions in agreement with ours, and also show that the results are in 
accord with the results of numerical simulations. 

The authors would like to acknowledge G. Kavoulakis and W. Ketterle for helpful discussions. 
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Captions of Figures 



FIG. 1. Frequency of the 0~-mode. This plot is obtained by taking the real part of the lower solution (y^FTojo) of the 
dispersion relation resulting from equating the A-dependent factor in Eq. (|l^) to zero in the limit A <C 1. The imaginary part 
is shown in Fig. 2. 



FIG. 2. Damping rate of the -mode, in the limit A ^ 1. 



FIG. 3. The hydrodynamic damping rates TJJ^ and in units of ojqt, as a function of A. The rate Fq. is associated with 
the 0^-mode while r!j+ corresponds to the 0'''-mode. 



FIG. 4. Damping rate versus mode frequency of the 0~-mode for A = \/8. The solid line represents the exact solution of Eq. 
( p^ and the dashed-dotted line is given by Eq. (|29|). The dotted line is the difference between the two upper curves (note the 
change of scale) . 



FIG. 5. Damping rate versus mode frequency of the -mode. The two curves are given by Eq. ( |29| ) for different parameters 
corresponding to the two experiments of Refs. [8,9]. The points indicate experimental results. 
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